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Abstract 

We synthesise traditional unstructured food webs, allometric body size scaling, trait-based mod- 
elling, and physiologically structured modelling to provide a novel and ecologically relevant tool for 
size-structured food webs. The framework allows food web models to include ontogenetic growth 
and life-history omnivory at the individual level by resolving the population structure of each 
species as a size-spectrum. Each species is characterised by the trait 'size at maturation', and all 
model parameters are made species independent through scaling with individual body size and size 
at maturation. Parameter values are determined from cross-species analysis of fish communities 
as life-history omnivory is widespread in aquatic systems, but may be reparameterised for other 
systems. An ensemble of food webs is generated and the resulting communities are analysed at four 
levels of organisation: community level, species level, trait level, and individual level. The model 
may be solved analytically by assuming that the community spectrum follows a power law. The 
analytical solution provides a baseline expectation of the results of complex food web simulations, 
and agrees well with the predictions of the full model on 1) biomass distribution as a function of 
individual size, 2) biomass distribution as a function of size at maturation, and 3) relation between 
predator-prey mass ratio of preferred and eaten food. The full model additionally predicts the 
diversity distribution as a function of size at maturation. 
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1. Introduction 

Food webs are typically modelled using unstructured species populations based on generalised 
Lotka-Voltcrra equations. This unstructured formulation ignores individual life-history by as- 
signing a fixed trophic position to all individuals within a species. In aquatic ecosystems this 
assumption is violated as fish offspring reside at a low trop hic level and grow during ontogeny 
through multiple trophic levels before reaching maturation ( Werner and Gilliaml . 1984 ). Along 



this journey, from the milligram range and up to several kilogram, fish change diet (as well as 
enemies) and conseque ntly exhibit life-history omnivory through preying on different trophic levels 
in different life-stages ( Pimm and Ricei . 1987 ). Thus the assignment of a unique trophic level and 



role (resource, consumer, predator, etc.) for species in unstructured models is incompatible with 
systems where ontogenetic growth and life-history omnivory are pronounced. In the cases where 



troph ic level of individuals within a species is positively correlated with body size (jJennings et al 



2002), individual size may be used as a proxy for trophic level. Models may therefore account for 
ontogenetic growth and life-history omnivory by resolving the size-structure within each species. 

A general framework for large food webs that includes the size-structure for all species must ful- 
fil a set of requirements. It should: 1) be generic in the sense that large species-specific parameter 
sets are not necessary, 2) be based on mechanistic physiological individual-level processes, where 
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paramet ers represent measurable b iological quantities, 3) resolve food dependent growth of indi- 
viduals ( Werner and Gilliairlfltil . 4) be practically solvable for species-rich systems over many 
generations, and 5) comply with empirical data on size-structured communities. In this work we 
develop a food web framework complying with these requirements by resolving the life-history of 
individuals within species by a continuous size-spectrum. We parameterise the model for aquatic 
systems as an example of a size-structured community with widespread life-history omnivory, but 
the framework may be parameterised for other system types (cf. Discussion). In fish communi- 
ties the most prominent empirical patter ns, which the model framework should comply with, are 
that i ndividuals exhi b it bip hasic growth ( Lester et al. . 20041 ). and the Sheldon community spec- 
trum. ISheldon et al.l ([19721 ) hypothesised that the community biomass spectrum, from bacteria 
to whales, as a function of body mass is close to constant. Empiric studies later showed t hat the 
biomass for fish indeed is close to constant or slightly declining as a function of body mass (jUrsinl . 



1982; Boudreau and 



Dickie, 1992) with the complication tha t heavily fished systems have a steeper 



Jennings et all . 120021 iDaan et all [2005) 



decline in biomass 

The importance of resolving ontogenetic growth and life- history omnivory has long been realised 
in fisheries science, wher e mechanistic individual-lev el size-structured food web models of fish com- 
munities were pioneered (|Andersen and Ursin , 1977). Independently, the physiologically struct ured 
population model (PSPM) framework ( Metz and Diekmann . 1986tlde Roos and Perssonl feoOl') has 
been developed in the field of ecology. While providing the ecological realism needed for a size- 
structured food web framework these approaches typically rely on large species-dependent param- 
eter sets, which must be reduced for the approaches to be useful as generic frameworks. 

Reduction to species-independent parameter sets has been achieved in unstructured mod- 
els of interacting popula tions by scaling of physiological and demographic rates with body size 
( Yodzis and Innesl . ll992h . By using body size as a trait t his approach has resu l ted in several sim- 
ple generic food web models for unstructured po pulations ( Loeuille and Loreaul 120051 : IVirgo et al 
2009: iBrose et all . l2006bl: iLewis and Lawl . l2007t ). 



In this work we combine the two approaches into one unified framework: We 1) use a physi- 
ological based description of individual life-history, and 2) use a single trait (size at maturation) 
to characterise each species while using trait and body size scaling to get one condensed species- 
independent parameter set. All processes are based on descriptions at the level of individuals, and 
interaction strengths among individuals are dynamic through the prescription of size-dependent 
food selection. This leads to a realised effective food web structure which depends on the emergent 
size-spectrum composition of all species. In this manner we synthesise a general framework that in 
a conceptually simple yet ecologically realistic way can be used to model food webs where the life 
cycle of individuals in each species is explicitly modelled from birth to reproduction and death. 

Our primary objective is the formulation and paramctrisation of the food web framework. 
Food webs generated by unstructured food web models may be analysed at the community level 
in terms of distributions of biomass across species and trophic levels. Trait-based size-structured 
food webs allow a more detailed analysis of the community level as well as enabling analysis on 
three additional levels of organisation: 1) at the community level, i.e., the distribution of total 
biomass as a function of body size of individuals regardless of their species identity, and the 
distribution of biomass and diversity as a function of the trait size at maturation, 2) at the species 
level, i.e., distribution of biomass as a function of size within a given species, 3) at the trait level, 
which in the case of a single trait equals the species level, and 4) at the individual level, i.e., 
distribution of size of food in the stomachs. Due to this added complexity of size-structured food 
webs, our secondary objective is to illustrate diversity and biomass distributions at different levels 
of organisation. To this end we generate an ensemble of food webs and analyse them in terms of 
distributions of average community size-spectra, species size-spectra, trait biomass distributions, 
and trait diversity distributions. Finally, we develop an analytical solution of the model framework, 
basically by assuming that the community spectrum follows a power-law (equilibrium size-spectrum 
theory, EQT). All distributions, except the diversity distribution, may be calculated from EQT, 
and we demonstrate general accordance between EQT and the results from the full food web 
simulations. The accordance between EQT and the food web simulations validates the simplifying 



assumptions behind EQT. EQT provides a "null-solution" to the size- and trait-distributions which 
may be used as a baseline expectation of the results of large size-structured food web simulations. 



[Figure 1 about here. 



2. Food web model 

The model is based on a description of the processes of food encounter, growth, reproduction, 
and mortality at the level of an individual with body mass to (Fig. [IJ. The model is based on 
two central assumptions: 1) Prey selection is determined at the individual level where individual 
predators select prey from the rule "big individuals eat smaller individuals" , a nd at the species level 
through introduction of species- s pecific size-independent couplin g strengths (jAndersen and Ursinl . 
Il977l : IWerner and Gillianl . ll984nEmmerson and Raffaelll 12004 ) . 2) In addition to species-specific 
coupling strengths, species identity is characterised by a single trait: size at maturation m*. Inter- 
actions among individuals are described by a food encounter process which leads to consumption 
by predators and mortality on their prey. Food consumption leads to growth in body mass, and 
when an individual reaches size at maturation it starts allocating energy for reproduction, as 
well as producing new offspring. Thus the model encapsulates the life-cycle of individuals from 
birth to maturity and death. 

Population dynamics of species i is obtaine d from individual growth q,;(m) and mo rtality (to) 



by solving the number conservation equation ( M cKendrickl Il926t Ivon Foersterl . [l959| ) 



dNj 
dt 



— (< 

dm V' 



N t = -fuNi. 



(1) 



The population structure of species i is described by the size-spectrum Ni(m,t), denoted Ni(m) 
to ease notation. The size-spectrum represents the volumetric abundance density distribution 
of individuals such that Ni (to) dm is the number of individuals per unit volume in the mass 
range [to; to + dm]. Similarly Bi(m) = mNi(m) denotes the biomass spectrum (biomass density 
distribution), and Bi(m) dm the biomass per unit volume in the range [to; to + dm]. The sum of 
all species' size-spectra plus a resource spectrum N R (m) is the community spectrum (Fig. [I}: 



N c {m) 



(2) 



The community spectrum represents the entire biotic environment providing individuals with food 
(from smaller individuals) as well as their predation risk from larger individuals. To include 
species-specific preferences each species i has its own experienced community spectrum: 



M(m) = e itR N R (m) + OijNjim), 



(3) 



where 8ij € [0; 1] is the coupling strength of species i to species j. Coupling strengths are 
independent of body size (cf. Discussion) since size-dependent food intake is described with a 
feeding kernel (below). 



2.1. Food consumption 

The consumption of food by an individual depends on the available food from the experienced 
community spectrum, on the volume searched per time, and on its functional response. The 
consumed food is assimilated and used to cover respiratory costs. Remaining available energy is 
used for somatic growth by immature individuals and for a combination of somatic growth and 
reproduction by mature individuals. 

We incorporate the rule of "big ones eat smaller ones" by assuming that predators have a 
preferred pr edator - prey mass r atio (PPMR). This assumption is inspired by stomach ana l yses o f 
marine fish (jUrsinl . Il973l . Il974h . and supported by stable isotope analyses (| Jennings et all 120011 ). 



The feeding kerne l descr i bing t he size preference for prey is prescribed with a normalised log-normal 
function (Fig.El lUrsinl . Il973l ): 



s(to p , m) = cxp 



In 



/(2a 2 ) 



(4) 



where m p is prey mass, to predator mass, f3 the preferred PPMR, and a the width of the function. 
The food available (mass per volume) for a predator of size to is: 



4>i(m) = / mpj\fi(m p )s(mp,m) drrip. 



(5) 



Encountered food (mass per time) is the available food multiplied by the volumetric search rate 
v{m) = 7m 9 , where o is a positi ve exponent signifying that larger individuals search a larger 



pc 

volume p er unit time (IWard. 1 19781). Satiation is described using a feeding level (number between 
and 1, Kitchell and Steward ." 1977 : Andersen and Ursinl . ll977h : 



v(m)(j)i (to) 
v(m)(f>i(m) + hm n ' 



(6) 



where hm n is the maximum food intake. Feeding level times hm n corresponds to a type II func- 
tional response. 

2.2. Somatic growth 

Ingested food f(m)hm n is assimilated with an efficiency a accounting for waste products and 
specific dynamic action. From the assimilated energy an individual has to pay the metabolic costs 
of standard metabolism and activity, kmP . Thus the energy available for growth and reproduction 
is: 

Ei(m) = afi(m)hm n - kmP . (7) 

Of the available energy a fraction -0(to, to*) is used for reproduction, and the rest for somatic 
growth: 

^1 — -ip(m,m^)jEi(m) Ei{m) > 
otherwise 



5,(771, to*) = 



(8) 



If the intake is insufficient to cover respiratory costs (£"; (to) < 0) growth is halted. Body size does 
not shrink when costs cannot be covered, instead starving individuals are exposed to a starvation 
mortality (see section 12. 4|) . The maximum asymptotic size M an individual can obtain is reached 
when all available energy is used for reproduction (t/>(M, to*) = 1). 

2.3. Reproduction 

In order to generate growth trajectories with biphasic growth the allocation rule ip(m, m*) has 
to change smoothly from around size at maturation to 1 at the theoretical maximum asymp- 
totic size M. The allocation rule ij)(m,m* ) is derived using two requirements: 1) that the size of 



gonads is proportional to individual mass (Blue weiss et all 119781). and 21 that s i ze at maturation 
is pr oportional to asymptotic size ( BevertonT 19921 : Froese and BinohlanL 200dl : iHe and Stewart! 
20011 ). To obtain an analytical solution as to how individuals allocate available energy to growth 
and reproduction we assume that the allocation rule is based on a constant feeding level /. Re- 
quiring allocation to reproduction to be proportional to individual mass, ip(m,m*)E(m) — k r m, 
gives 4'(' m i rn *) = k r m/E(m), where E(m) = afhm n — km p denotes the available energy when 
feeding level is constant. The factor k r is found by the second requirement through ij){M, to*) = 1: 
k r = E{M)/M where M = m*/?7*. The allocation can thus be described as: 



where the term in the square brackets is a smooth step function switching from zero to one around 
the size at maturation (u determines transition width). 

The exponents of maximum consumption and standard metabolism are close to equal (cf. | Appendix E| 
and Discussion). In the limit of n = p the available energy for growth and reproduction becomes 
E(m) = hm n where H = afh — k. This gives: ip(m,m*) = [1 + {m/ra st )~ u ]~ 1 (rj st ra/m t ) 1 ^ n , 
meaning that the juvenile growth pattern is g = hm n whereas adults grow according to g = 
Hm n — h(m* J l i]*) n ~ 1 m. Thus the growth m odel is a biph a sic gr owth model where adults follow von 



Bertalanffy growth curves as advocated bv lLester et al. (2004). 

The total flux of offspring is found by integrating the energy allocated to reproduction tj)(m, m*)Ei(m) 
over all individual sizes: 

Ri = / Ni{m)^(m,m*)Ei{rn) dm, (10) 

2m J 



where mo is the egg size, e the efficiency of offspring production ( | Appendix C I , and 1/2 takes into 
account that only females spawn (assuming equal sex distribution). Reproduction determines the 
lower boundary condition of ([1]) for the size-spectrum of the species: 



i(mo,m*)Ni(mo) = Ri 



(11) 



24. Mortality 

The mortality rate n(m) of an individual has three sources: predation mortality /j, p (m), star- 
vation mortality /x s (m), and a small constant background mortality /Zb(m*). The background 
mortality is needed to ensure that the largest individuals in the community also experience mor- 
tality as they are not predated upon by any individuals from the community spectrum. 

Predation mortality is calculated such that all that is eaten translates into predation mortalities 
on the ingested prey individuals ( | Appendix A ): 



P P A m p) = / s ( m p> m )(l _ fj(rn))v(m)6j ii Nj(m) dm. 



(12) 



When food supply does not cover metabolic requirements kmP starvation mortality kicks in. 
Starvation mortality is proportional to the energy deficiency km p — af(m)hm n , and inversely 
proportional to lipid reserves, which are assumed proportional to body mass: 



/i a (m) 



Ei(m) > 

■^'( m ) otherwise 



(13) 



Mortality from other sources than predation and starvatio n is assumed constant within a species 
and inversely proportional to generation time ( Petersl 1983 ): 



n — 1 

lib = Mom, • 



(14) 



2.5. Resource spectrum 

The resource spectrum Nn(m) represents food items which are needed for the smallest individ- 
uals (smaller than /3mo). The dynamics of each size group in the resource spectrum is described 
using semi-chemostatic growth: 



dN R (m,t) 
dt 



r^m H 



nm - Nn(m,t) - n P {m)NR{m,t), 



(15) 



where rgm^ 1 is the population regeneration rate ( Fenchell . 1974 : Savage et al. . 2004 ) and nm~ x 
the carrying capacity. We prefer semi-chemostatic to logistic growth since planktonic resources 
rebuild from depletion locally due to both population growth and invasions. 



[Table 1 about here.] 



2.6. Derivation of parameters 

Each species is characterised by a single trait, size at maturation to*, and a species-independent 
parameter set is achieved through scaling with body size to and to*. The model is parameterised 
for marine ecosystems using cross-species analyses of fish communities (|Appendix E and Table [1]). 



The constant 7 for the volumetric search rate is difficult to assess (Appendix E|. However, 
since the feeding level /(to) of small individuals is determined solely by the amount of encountered 
food from the resource spectrum, we may use initial feeding level fo as a physiological measure of 
food encounter; fo is defined as the feeding level resulting from a resource spectrum at carrying 
capacity. The initial feeding level is used as a control parameter for food availability (enrichment), 
through which the value of 7 can be calculated ( | Appendix D ) : 



7i(/o) = n f \ FT A ' (16) 

where it is noted that 7 will be species dependent if species have different coupling strengths to 
the resource. 

A critical feeding level f c can be formulated as the feeding level where all assimilated food is 
used for metabolic costs (using values from Tabled]): 

fc = \rnP~ n = A « 0.2. (17) 
an an 

Individuals can only grow and reproduce if / > f c . Assuming that individuals experience an 
average feeding level /, the growth (JSj) of juveniles is g = hm n (for n — p). The parameter 
h = ahf — k can be estimated through the rela tion between observed von Bertalanffy growth rate 
and asymptotic size yielding h w 10 g 25 /year ( Andersen et all 120081) . This allows an estimation 



of the expected average feeding level of individuals in the field (Table [T]): 

7=^-0.4, (18) 
an 

i.e. around twice the critical feeding level. As the initial feeding level fo is calculated from a 
resource spectrum at carrying capacity, the realised feeding level in the model will be smaller than 
fo. A value of f = 0.6 was seen to give realised feeding levels around 0.4. 



3. Methods 

Stable food webs are constructed using the full dynamic food web model with random coupling 
strengths 9ij. For each run, 30 species are assigned with m* evenly distributed on a logarithmic size 
axis (to* € [0.25 g; 20kg]), random 9ij matrices (mean 0.5), and a common 9i t R = 0.5 coupling to 
the resource spectrum. Numerical integration is performed by standard finite difference techniques 
( |Appendix G[ ). Food webs are simulated in 10 consecutive intervals covering 300 years each, where 
species with a biomass less than 10 _30 g/m 3 are eliminated after each interval. To eliminate food 
webs that still have not reached the final state each community is integrated for additional 500 
years and discarded if any species has an absolute population growth rate larger than 1 logarithmic 
decade per 100 years. To ensure that each food web in the final ensemble spans multiple trophic 
levels we only retain food webs where at least one species has m* larger than 2.5 kg. For statistics 
we use the mean of the last 250 years of the simulation with time steps saved in 0.1 year increments. 
In this manner 204 food webs having a total number of 1016 species were collected. Each web 
contained between 2 and 9 species with a mean of 4.98 species. 

We analyse the generated food webs in terms of distributions of average community size- 
spectrum, species size-spectra, trait biomass distributions, and trait diversity distributions. Ad- 
ditionally we demonstrate the importance of distinguishing between what an individual prefers to 
eat and what is actually ingested (i.e. found in its stomach) by showing how emerging PPMRs 
vary with food availability and differ from preferred PPMRs. 



An approximate steady-state solution to the food web model which neglects the dynamics of 
reproduction can be found using two assumptions: 1) all species consume food and experience 
mortality from a scaling community size-spectrum N c = K c m _A , and 2) constant feeding level /, 
which implies equal species coupling strengths O^j — 8. Whereas the food webs in the full model 
are based on a discrete set of to* , the analytical solution considers to* as a continuous distribution. 
The procedure for deriving the analytical sol ution is similar to the derivation of equilibrium size- 
spectrum theory (jAndersen and Beverl 120061 ) . but the results are slightly different as standard 
metabolism is taken explicitly into account here. The food encountered by an individual is found 
using assumption 1): v(m)<f)(m) = 7m 9 J N c s(mp,m)m p dm p oc m 2 ~ x+q . The feeding level is 
calculated from and the requirement that it is constant (assumption 2) leads to a constraint 
on the exponent of the community spectrum: A = 2 + q — n. Feeding with a constant feeding level 
generates a predation mortality of fi p = a p m n ~ 1 ( | Appendix A ) . The size-spectrum of juvenile 



individuals is found as the steady state solution of (JTJ) using the above predation mortality and 

g = hm n fcf . (IF. 11) 1: N (rn,m*) = K(m*)m~ n ~ a , w here a = a p /h is the physiological level 

of predation (jBeveri Il989t lAndersen and Bever . l2006h . which can be calculated as a ss //(/ 



Qpin-q-Xja _ g6 ( [Appendix Bj ). The constant «(m») is found from the requirement that the 
sum of all species spectra should equal the community spectrum. Assuming a continuum of species 
the requirement can be written as J N(m, to*) dm* = N c (m) which leads to k(to*) oc m* n q 3+a 
(Fig. [1} . This approximate solution of the model will be referred to as equilibrium size-spectrum 
theory (EQT), and it will be compared to the output of the complete dynamic food web model. 

In dynamic models, as in nature, the lifetime reproductive success (fitness) has to be Ro = 1 
for all coexisting species. Since EQT does not consider the boundary condition (TTTI) life- time 
repro ductive success becomes a function of size at maturation: R oc m\" a (jAndersen et all 
120081) . One solution to making R independent of to* is to set a = 1, but that breaks the above 
employed mass balance between growth and mortality used to calculate a. Due to the Rq ^ 1 
inconsistency in EQT we have a specific focus on the realised values of a when comparing food 
web simulations with EQT predictions. To examine how the regulation of Rq occurs in the full 
food web model Ro is split into two factors: 1) the probability of surviving to become adult, and 
2) lifetime reproduction per adult ( | Appendix F[ ) : 

_ N(m)g(m,m*) 
P™^- N{mo )g(m 0i m*y {U) 

u ( \ f M ip(m,m*)E(m) 

-ttadult(m*) = / Pm.^rn ; s dlJl. (20) 



g(m, m*) 

Survival probabilities and reproductive outputs in the food web simulations are compared 
with EQT predictions, which are calculated by inserting the EQT size-spectra into (fT9|) and ([20]) . 
Juvenile growth is g oc to™, which gives y mQ -+m, oc m~ n ~ a m™ = m~ a and i? a duit oc m». 

[Figure 2 about here.] 



4. Model predictions 

4-1. Growth trajectories 

In unstructured models fluctuations are manifested as oscillations in the biomass of species, 
whereas the oscillations in structured models stem from oscillations in the size-spectrum compo- 
sition. Such oscillations give rise to fluctuating feeding levels as individuals encounter different 
levels of food in different life-stages (Fig. (5]a). Variations in feeding levels between species and as 
a function of individual size lead to different emergent growth trajectories (Fig. (21b). The growth 
trajectories roughly follow the biphasic growth curve that is obtained if the feeding level is assumed 
to be constant. 



[Figure 3 about here.] 
[Figure 4 about here.] 



4-2. Biomass structure 

By pooling species from each food web into logarithmic evenly distributed to* groups, and 
summing the size-spectra in each group, a size-spectrum is obtained for each to* group. Next, 
the logarithmic average of m* groups across all food webs is performed to produce the average 
size-spectra of a m* group (Fig. [3|). Average community biomass spectrum N c (m)m follows the 
EQT prediction of a slope of 1 + q — n = 1.05, meaning that the biomass in logarithmically 
evenly sized size-groups, f N c (m)mdm : is a slightly declining function of body mass. The 
community spectrum oscillates around the EQT prediction due to a trophic cascade initiated 
by a superabundance of the la rgest predators which do not experience any predation mortality 
( Andersen and Pedersen , 2010f ). The peaks of the oscillating pattern are roughly spaced by the 



preferred PPMR. Biomass density within species is constant until individuals reach the end of the 
resource spectrum, and larger individuals, > 0.1 g, have a biomass spectrum slope flatter than 
that of the community spectrum (Fig. [3]). Thus, in contrast to EQT, the dynamic model produces 
species size-spectra that cannot be described as power laws. The number of small individuals 
is inversely related to size at maturation. The scaling of offspring abundance can be calculated 
using EQT as Nq oc J^™* K{m'^)dm'^ oc mt n ~ q 2+a , which fits the simulated results well for a = 1 
(Fig. El inset). 

The distribution of species biomass as a function to* can be calculated from EQT as: 

/■cm* pM 

-6(771*)= / / N(m, m'^mdm dm* oc m* ~ q . (21) 

•J m* J itiq 

As 77 and q are almost equal the biomass distribution -B(to*) as a function of to* is almost constant. 
This result is also borne out by the dynamical simulations (Fig. |4ja) with some variation due to 
uneven species distribution along the to* axis: peaks occur in species diversity separated by the 
preferred PPMR /3 (Fig. Hb). 

[Figure 5 about here.] 

4.3. PPMR and feeding level 

The realised mean PPMR can be derived when prey concentrations are known: Af(m p )s(m pi to) 
is the prey size distribution encountered by a to sized predator. Mean prey size encountered by a 

Jq 00 m p j\f(m p )s(m pi m)dm. 
J °° N(m p )s(m p ,m)t 

size to divided by the mean prey size 



to sized predator is ° roo p . r , , , — ^Vf — - . The realised mean PPMR is calculated as the predator 

' j A/ (m p )s(m p ,m)dm p 1 



PPMR(m) = ^^Mrn P ^)dra r 
J rn p Aj{rn p )s(rn p ,m)dm p 

Realised mean PPMR is always larger than the preferred PPMR j3, due to higher abundance of 
smaller prey items (Fig. [S]a). The realised mean PPMR calculated from EQT (using Af oc m~ A ) 
is exp[(A — 3/2)a 2 ]/3 1.7/3. Realised PPMR from the simulations oscillate around this value due 
to the fluctuations in the community spectrum (Fig. [3J). 

As individuals grow to a size larger than /3too they switch from eating food in the resource 
spectrum to feeding on other species. This leads to a decrease in the feeding level from /q = 0.6 
to about 0.45. The oscillations in feeding level increase in magnitude as body size increases due 
to larger fluctuations in prey availability (Fig. [5]b). Many large individuals periodically have a 
feeding level below the critical feeding level f c (where starvation kicks in) since prey items in the 
preferred size range become scarce, which results in ingestion of smaller food items and therefore 
increasing PPMR. 



[Figure 6 about here.] 



4-4- Reproduction and survival 

Lifetime adult reproduction calculated from the simulated food webs fit the EQT prediction 
since it scales linearly with m* (Fig.[6]b). The probability of surviving to a given size is independent 
of m*, as the survival curves of the different m* groups lie on top of one another (Fig.[6]a). Survival 
to m* scales inversely with m* (i.e. a = 1, crosses in Fig. Ela), which ensures that Ro is constant. 
However, if the a = 1 scaling of survival to m* is extrapolated to mo it is seen that it does not 
intersect Pm ->m = 1- Instead the survival curves change slope between mo and around 10 _1 g 
where predation mortality starts to dominate due to an abundance of fish individuals in the same 
order of magnitude as the resource spectrum, which is intensified by reduced growth stemming 
from food competition (Fig. [3]). In summary survival does not scale with m~ a as predicted by 
EQT. Instead adult survival scales with m^ 1 (i.e. a = 1) whereas individuals smaller than rj 0.1 g 
have a higher survival (i.e. a smaller scaling exponent). 



5. Discussion 



We have developed a generic food web framework suitable for analysing systems of interacting 
size-structured populations. The framework increases ecological realism compared to traditional 
unstructured food web models by explicitly resolving the whole life-history of individuals, but 
maintains simplicity by describing species with only one trait: maturation size m». Remaining 
parameters are made species-independent through inter- and intraspecies scaling with m* and 
body mass m. The productivity of the system is characterised by one parameter, the initial 
feeding level fo- Feeding behaviour is assumed to be determined by a feeding kernel with a fixed 
preferred PPMR (big individuals eat small individuals), multiplied by a species-specific coupling 
strength. 

Only characterising the life-history and feeding preference of individuals of a species by body 
mass m and trait m* is clearly a simplification, but contemporary knowledge sugges ts that a large 



mass m ana trait m* is clearly a simplification, out contemporary knowledge suggest s tnat a large 
part of the individual bioenergetics related to growth ( Petersi 19831) and reproduction ( Blueweiss et al.l 
19781 ) indeed can b e described by such scaling. Additionally it is well-known that predators often 



outsize their prey (jBrose et al 
small ones" . 



2006al ) which justifies the use of the generalisation "big ones eat 



5.1. Model architecture 

The model was parameterised from cross-species analyses of fish communities, since aquatic sys- 
tems constitute a group of strongly size-structured ecosystems. Other less strongly size-structured 
taxa can be modelled as well through reparametrisation and by allowing each species to have 
its own offspring size mo.i- Additionally, the description of how available energy energy is divided 
between growth and reproduction may have to be reformulated since animals in other taxa may ex- 
hibit determinate growth. Determinate growth can be modelled simply by replacing the allocation 
function ^ with only the part within the square brackets. 



T he proposed modelling framework is similar to phys i ologic ally structured models (jAndersen and Ursin 
Metz and Diekmannl . Il986t Ide Roos and PerssonlEool . and as these based on mechanistic 



1977 



individual-level processes. Our contribution is to employ a trait-based description of species iden- 
tity, and a formulation of food preference which is split into a size- and a species-based contribu- 
tion, which renders the developed framework useful as a generic food web framework. Recently 
the PSP M approach has been reduced to a stage-structured model which facilitates multi-species 



studies (|de Roos et all l2008aj ): however this is achieved at the cost of collapsing continuous size- 



structure to a discrete stage-structure. A first step towards multi-species PSPMs was carried out 
with an intra-guild predation model, which sho wed that obtaining species co existence between 



two size-structured populations is a difficult task (jvan de Wolfshaar et all 12006); a result which is 



probably due to insufficient ecological differentiation of the two species. In the proposed framework 
the trait maturation size provides a simple and logical way of representing ecological differentiation 
of species, whereas this differentiation in PSPMs is less clear due to large species-specific parameter 



sets. Additional ecological differentiation and heterogeneity are obtained by also including food 
web structure in the form of species coupling strengths. 

An alternative approach to model a size-structured com munity is the community size-spectrum 

models (|Silvert and Platu 11980c iBenoit and Rochetl . 12004) . In these models the comm unity is rep- 

resen ted by a community size-spectrum of all individuals irrespective of species identity (jSheldon and Parsons , 

)67j ). As with the physiologically structured models these are based on individual- level descrip- 
tions of life-history. The community spectrum approach has the drawback that species are not 
resolved, as all individuals are lumped together into one s pectrum. Their advanta ge is their ability 
to make community-wide predictions with si mple means (iBlanchard et all [2 009) simil arly to the 
mean- field theory in unstructured food webs (jMcKane et all l2000t IWilson et all 120031 ). 

A central element in the model is the division of energy between somatic growth and repro- 
duction through the allocation function ip(m,m*). As in PSPMs our bioenergetic model is a net- 
production model where it is assumed that metabolic costs are covered with highest priority after 
which the remaining energy can be used for grow th and reproduction. PSPMs are formulated either 
with one state variable: individual body weight ( Kooiiman and Metzl .ll984; Cla essen and de Roos 



2003), or with two state variables: somatic weight and reserve weight (|de Roos and Perssonl 120011 ). 



In the latter case energy is divided between the two states such that the ratio between the two 
state variables is aimed to be constant, and accumulated reserves are used for reproduction at 
the end of the growing season. In the case with only one state variable surplus energy is divided 
between somatic growth and reproduction with a fixed ratio (re-rule). When using the re-rule the 
maximum asymptotic size any species individual can obtain is M+ where intake ahf(M + )M™ 
equals the metabolic costs kM+ - meaning that all species would obtain the same asymptotic size 
if parameters are species independent as in our formulation. M+ is very sensitive to the precise 
valu es of n and p, and they can therefore only be regarded as poor detcrminators for asymptotic 
size (j Andersen et all 12008). Our model deviates from the single-state PSPMs in this partitioning 
of energy, as we assu me that mature i ndivid uals allocate an amount proportional to their body 
size for reprod uction dBlueweiss et al. . 1978), and that asymptotic size depends o n the trait size 
at maturation (jBevertonl . Il992t iFroese and Binohlanl . 2000; iHe and Steward . l200lh . This ensures 
that the ratio between gonad size and somatic weight is constant within a species, which is sim- 
ilar to the partitioning rule used in two-state PSPMs. The allocation function is derived under 
the assumption of a constant feeding level throughout adult life. Even though the feeding level 
is assumed constant, the actual allocation still vary depending on the actual food availability, as 
tp(m,m*) only determines the fraction of available energy allocated to reproduction. An alterna- 
tive way to derive ip(m, to*) is to let it depend on actual feeding levels. This assumption, however, 
would imply that individuals adjust their allocation to reproduction such that asymptotic size 
is always reached. This does not seem plausible as individuals in resource scarce environments 
probably obtain smaller maximum sizes, and therefore we find the most reasonable assumption to 
be that of a constant feeding level. The exponents n and p are close to equal in nature, and for 
n = p the ener gy allocation f unctio n leads to biphasic growth where adults follow von Bertalanffy 
growth curves f Lester et all 12004 ) • We fixed the yearly mass-specific allocation to reproduction 
(yearly gonado-somatic index, GSI) to be independent of individual body size within a species. The 
obtained form of ip(m, m*), how ever, yields a to"" 1 scaling of yearly GSI across species, which is 
consistent with empiric evidence (| Gundersonl . 1 1 99 71) . This means that the form of ip(m, to*) implies 
a trade-off between to* and the mass-specific reproduction: large to* species can e scape predation 



mort ality via growth by paying the price of a lower mass-specific reproduction (jCharnov et al 



200lh . When the exponents n and p differ, growth will still be biphasic and adul t growth curves 



will be similar to von Bertalanffy curves (see also Andersen and Pedersen . 2010f) . In conclusion 
the derived allocation rule leads to realistic growth patterns. 



5.2. Food web structure 

Food web structure is the most essential part of a food web model, and in principle two ap- 
proaches can be taken to obtain a structure for a dynamic food web model: a top-down and a 
bottom-up approach. 



The top-down method generates food web matrices from the desire d number of species and con - 
nectance using a static mode l (stochastic phenom e nological mode ls: Cohen and Newmanl ( 1985^ : 
Williams and Martinez! (|2000l ): fCattin et al.l d2004h : lAllesina et all (|2008l ). or more me chanistic ap- 
proac hes involving phylogenetic correlations (jRossberg et all 120061) or foraging theory ([Petchev et al 
20081 )). Next, the food web matrix is used to drive a dynamic model, which is simulated forward 
in time where some of the initial species will go extinct, and the remaining species set can be used 
for analysis. Note that in addition to a decreased species richness in the final community o ther 
food web statistics as e.g. the final connectance may differ as well ( Uchida and DrosselL 2007 ). 

In the bottom-up approach link strengths are determined from ecological relations, such as 
e.g. a predator-prey feeding kernel: if the prey fits into a certain size range relative to the preda- 
tor siz e, then interaction occurs between the nodes with a strength dete rmined by the feeding 



kernel ([Loeuille and LoreauL 12005c I Virgo et all 12006c iLewis and Lawi . 12007) . Predator preferences 



depend, in addition to ecological characters, on evolutionary hist ory and recent approaches add 
this component of phylogenetic correlations (jRossberg et al . 120081 ). 



For size-structured food webs a top-down algorithm for generating realistic food web matrices 
does not exist. This is due to lack of data describing the three dimensional interaction matrix 
- dimension one and two is respectively predator and prey identity as in the classic interaction 
matrix, and the third dimension is predator/prey body size. Thus one is confined to the bottom- up 
approach and/or random interaction matrices. In this study we use the bottom-up approach to 
prescribe interac t ions t o obey the pattern of "big ones eat smaller ones". Life-history omnivory 
( Pimm and Ricd . Il987l) is therefore naturally incorporated in size-structured food webs through 
the use of a feeding kernel. To obtain an ensemble of different communities we use the top- 
down approach of a classical two-dimensional predator-prey interaction matrix - that is we assume 
that regardless of size individuals within a species have equal potential maximum link strength 
(coupling strengths in our model) to another species. As no top-down method exists for generating 
this matrix we use random matrices. The actual link strength is the product of the coupling 
strength and the feeding kernel, meaning that link strengths indeed are dynamic as they depend 
on the size-structure of both prey and predator. 

As we generate food webs from a fixed initial pool of only 30 species and use a random matrix 
as coupling matrix we only obtain small food webs (maximum: 9 species). However, it should be 
noted that the number of resource species the resource spectrum represents is not included. To 
obtain larger foo d webs a larger species pool is needed along with a sequential assembly algorithm 
(jPost and Pimml . Il983l ). and a better method for obtaining coupling strengths between species. 
Our primary interest in the food web analyses has been the size- and trait-structure of food webs 
with a finite number of species, and how these compare with EQT predictions, which are based 
on the premise of a continuum of species. The general correspondence with EQT indicates that 
the broad-scale patterns are relatively insensitive to how the species-specific coupling strengths 
(i.e. food web structure) are specified. Still, an interesting follow-up study would be focused on 
the coupling matrix structure, which may more generally be size-dependent, and how the effective 
food web structure that emerges from the coupling strengths and feeding kernel compares with 
empiric food webs. 



5.3. Community structure 

We generated an ensemble of size-structured food webs and used averages over these to make 
general predictions of the structure offish communities, in particular the size-structure of individual 
populations, and how these spectra "stack" to form the community size-spectrum. In accordance 
with EQT we find the community spectrum to scale with A = 2 + q — n«2 meaning that the 
distribution of biomass as a function of individual body size is close to constant when individuals are 
sorted into logarithmically evenly sized bins. This prediction means that the biomass of individuals 
between e.g. lg and 10 g is the same as those present with bo dy sizes between 1kg and 10 kg, in 



accordance with the Sheldon hypothesis ( Sh eldon et al. ■ Il972h 



The distribution of biomass as a function of m* is predicted to be almost independent of 
m* in accordance with EQT. The result is reminiscent of the Sheldon hypothesis, and it can be 



formulated as an extension of the Sheldon hypothesis: "The total biomass of individuals ordered 
in logarithmically spaced groups of their maturation size is approximately constant". This means 
that the total biomass of all species with to* between 1 and 10 g is approximately the same as 
that of species with to* within 1 to 10 kg. This prediction is a novel extension and could be tested 
by size-based field data. In contrast to EQT the dynamic framework also provides predictions on 
the distribution of the number of species as a function of to* . Species tend to cluster in groups on 
the to* axis separated by a distance corresponding to the preferred PPMR j3. This clustering is 
partly a reflection of the use of a fixed value of /?; more diversity in feeding strategies (i.e. different 
/3) would probably smoothen the species distribution as well as making the feeding level more 
constant. 

The size-spectra of individual species do not to follow power laws as predicted by EQT since 
there is a change in spectrum slopes from small to medium sized individuals. This difference stems 
from different scaling relationships for the survival probability of small and larger individuals. The 
less steep slope in survival for small individuals is due to a proportionally low mortality rate caused 
by their low abundance relative to similarly sized resource items. Incorporating mortality from 
the resource spectrum on the smallest individuals may thus result in a single survival probability 
scaling. The probability of surviving to to* scale as m~ a for a physiological predation constant 
value of a = 1, which is conflicting with the value a = 0.86 predicted by EQT. The discrepancy 
about the value of a highlights an inconsistency within EQT: Enforcement of mass-balance between 
growth and predation leads to a = 0.86, while the reproductive boundary condition can only be 
fulfilled if a = 1. The full food web simulations demonstrate that both the scaling of surviving to 
to* and the scaling of the number of offspring are best predicted by a value of a = 1 . This indicates 
that when EQT predictions depend on a, the value a = 1 should be used even though that breaks 
mass conservation in EQT. 

Lastly we demonstrate that realised PPMRs (i.e. PPMRs based on ingested prey) emerge in 
the model. Average realised PPMR is always larger than the preferred PPMR /3 since smaller prey 
items are more abundant than larger ones. It is found that the realised PPMR is proportional 
to the preferred ratio (PPMR = 1.7/3). Model predictions show that realised PPMR oscillates 
around this value due to fluctuations in the average community spectrum. PPMR displays large 
fluctuations with size demonstrating that determination of PPMR from single measurements is 
problematic due to high prey abundance sen sitivity. Empirical findings show that realised PPMRs 
increase with body size (| Barnes et al.l . |2010|) , but one should be careful about concluding that the 
preferred PPMR (which we put into models) shares this size scaling, since relative abundances 
may cause the increase rather than actual behavioural prey preferences: even though we have a 
fixed preferred PPMR our model predicts that realised PPMR is an increasing function of body 
size. 



5.4- Conclusion and outlook 

The proposed food web framework increases ecological realism in food web models as it resolves 
the complete life- history of individuals by representing the size-structure of each species with a size- 
spectrum. More specifically the framework complies with five requirements of (cf. Introduction): 
1) being generic with few parameters, 2) being mechanistic and utilising individual-level processes, 
3) including food dependent growth, 4) being practically solvable for species-rich communities, and 
5) complying with data on community structure and individual growth curves. 

Trait-based size-structured food webs can be examined at four levels of organisation: at com- 
munity level, at species level, at trait level, and at the individual level. We generated empirically 
testable hypotheses of mainly biomass distributions at different levels of organisation. 

By assuming a power law community sp ectrum and a constant feed ing level the full dynamic 
model can be simplified to an EQT model ( Andersen and Beyerl 120061 ). Correspondence of pre- 
dictions by EQT and the full model justifies the use of the simplifying assumptions. EQT is a 
powerful analytical tool that in a simple manner yields insight to e.g. the biomass distributions 
within and across species in size-structured food webs. However, as EQT assumes steady-state, 



the study of emerging effects, e.g. diversity and responses to perturbations, has to be conducted 
with the full model. 

The PSPM framework has showed existence of alter native stable states where single populations 



can exist with diffe r ent siz e-structure compositions (|de Roos and Perssonl l2002t iPersson et al 
" 2008blh It is an open question whether such alternative states become more 
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widespread or if they disappear when more species interact with each other. This question is 
important since it tells whether such alternative states are expected to occur frequently or rarely 
in nature, and consequently whether exploitation can easily induce shifts between states. An 
important future challenge is thus to study the possibilities of multiple states in complex food 
webs - not only of single individual populations, but of the ecosystem as a whole. The proposed 
framework allows exactly this kind of studies since it provides a full ecologically realistic but 
conceptually simple model of size-structured ecosystems. 

Natural future extensions of the model could be to allow the species coupling strengths to 
be size-depend ent and make coupling strengths depend on vulnerability and forageability of prey 



and predators (jRossberg et all 120081 ) as well as on the spatial overlaps of the interacting species. 



Adding this extra level of mechanistic realism would allow the framework to be useful for studying 
ecosystem consequences of spatial changes of species populations, which could be driven by climatic 
changes. 
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Figure 1: Illustration of the community model: resource spectrum (with cut-off m cu t) and four species size-spectra 
(m* = 2.5, 25, 250, and 2 500 g) having off-spring of size m.n. Shaded regions mark the spawning stock from 
maturation size m* to maximum asymptotic size M. The sum of all spectra gives the community spectrum. The 
spectra shown are the steady-state solutions from Andersen and Beyer ( 2006 ) . Inset shows how individuals feed on 
smaller prey using a feeding kernel with a preferred PPMR of /3. 
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Figure 2: Example of an eight species cyclic state, a) Top: Feeding levels of the species along with the min/max 
(light grey) and the 25%/75% (dark grey) percentile values of the time-series. Dashed lines indicate initial /n and 
critical f c feeding level. Bottom: The time-average of the resource and species spectra along with the community 
spectra (thick dashed). The idealised community spectrum re c m -A (thin dashed) and the species maturation sizes 
m- t (circles), b) Time averaged growth curves for the species (thin lines) along with the biphasic growth curve J8jl 
for a fixed feeding level that equals 75% of the time and size averaged feeding level experienced by the species (thick 
line). Growth curves are normalised with realised asymptotic size (j/-axis) and generation time (x-axis) to enable 
comparison. 
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Figure 3: Mean species biomass spectra (grey lines) when species are divided into 5 logarithmic evenly distributed 
m* groups. Also shown is the total mean community biomass spectrum (thick line), and the EQT community 
biomass spectrum re c m 1— A (dashed). Inset shows how offspring abundance (No) scales with m* (data pooled in 25 
log groups). Expected EQT scaling of Nq is shown for a = 0.86 (dashed) and a = 1 (solid). 
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Figure 4: a) Distribution of biomass in different m* groups. The expected distribution B(m t ) oc m™~ 9 is illustrated 
with the solid line. Dashed line shows biomass per species (bar values divided with bar values in b)). b) Mean no 
of species as a function of m* . Species are pooled into 16 logarithmic groups. 




Figure 5: a) Realised PPMR from each food web (grey), mean realised PPMR across all simulations (thick black), 
realised PPMR prediction from equilibrium theory (thin black), and preferred PPMR /3 (thin dashed), b) Feeding 
level from each species (grey), and mean feeding level (thick black). 
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Figure 6: The components of the expected lifetime reproductive output, a) Probability of surviving to maturation 
size m* (crosses) along with survival curves throughout life (grey line) for data pooled into 7 logarithmic groups. 
Expected p <x scaling from EQT is shown for a = 0.86 (dashed) and a = 1 (solid), b) Lifetime adult 

reproduction for all species (crosses) along with EQT -R a dult <x m * expectation (line). 



T able 1: Default parameter values for a temperature of 10° C j Appendix E| . 
Symbol Value Units Parameter 
Individual growth 



fo 


0.6 


- 




Initial feeding level 


a 


0.6 


- 




Assimilation efficiency 


h 


85 




'"/year 


Maximum food intake 


n 


0.75 


- 




Exponent for max. food intake 


k 


10 


g 1 " 


_p /year 


Std. metabolism and activity 


P 


0.75 


- 




Exponent of std. metabolism 


13 


100 


- 




Preferred PPMR 


a 


1 


- 




Width of feeding kernel 


q 


0.8 


- 




Exponent for search volume 


Reproduction 








m 


0.5 


mg 




Offspring mass 


V* 


0.25 






m* rel. to asymptotic mass M 


e 


0.1 






Efficiency of offspring production 


u 


10 






Width of maturation transition 


Mortality 












0.1 






Fraction of energy reserves 


Mo 


0.84 


g 1 " 


""/year 


Background mortality 


Resource 


spectrum 








K 


5- 10~ 3 


g A 


-Vm 3 


Magnitude of resource spectrum 


A 


2 - n + q 






Slope of resource spectrum 


TO 


4 


g 1 " 


_p /year 


Regeneration rate of resources 


m cut 


0.5 


g 




Upper limit of resource spectrum 



Appendix should be available as online material. 



Appendix A. Derivation of predation mortality 

Predators with a size between m and m + dm have a food intake rate of s(m p , m)f(m)hm n 8N(m)dm 
for m p sized prey. The total density of food available from all prey sizes to the predators in [m;m + dm] 
is (j)(m) (0), meaning that the mortality experienced by a m v sized individual is: 

= W s{m p ,m)f 3 (m)hm^N 3 jm) ^ ^ 
J <Pj (m) 

The maximum food intake may be expressed as a function of f(m), v(m), and <f>(m) via such that 
the predation mortality can be written as: 

Mp,i( m p) = / s(m p ,ro)(l — fj(m))v(m)9j t iNj(rri) dm. (A.2) 
i 

By using the EQT assumptions of constant feeding level and a power law community spectrum (cf. 

section [3]) the mortality reduces to /x p (m p ) = 6(1 — /) J s(m p ,m)v(m)K c m~ x dm, which can be solved 
analytically: 

Mp( m p) = "pnip" 1 ' (A-3) 
where a p = 6(1 - J)V2^K c ja(3 1+q - x exp [\<r 2 (l + q- A) 2 ] . 



Appendix B. Available food and the physiological level of predation a 

Using the EQT assumption of a power law community spectrum allows calculation of the available 
food density <j)(m) — 8 J s(m p ,m)K c m p x m p dm p : 

<j>(m) = a^9n c m 2 ~ x , (B.l) 

where a 4, = V2^af3 x ~ 2 exp [\o 2 (2 - A) 2 ] . 

Using the EQT assumption of constant feeding level yielding A = 2 + q — n allows us to write 9k c = 
fh/(a^(l — /)) by rearranging the expression of the feeding level ((6]). Using this and the definition 
of h allows writing a p = c(h + fc)/? 2 " -9-1 /a where c = exp [±cr 2 ((1 + q - A) 2 - (2 - A) 2 )] = 1.03 w 1. 
Ultimately using the definition of f c allows writing the physiological level of predation a = a p /h as: 

a = c^—^-^la. (B.2) 



Appendix C. Calculating efficiency e of offspring production 

The efficiency of turning energy into offspring is denoted e. It includes losses due to behavioural aspects, 
pre-hatching mortality, and that the energy contents in gonadic tissue is higher than in somatic tissue. It 
is a quantity that is difficult to measure, but for n = p its value can be derived. 

The energy (in units of mass) routed into reproduction (for n = p) is ip(m, m,)hm n where h = afh — k. 
The energy of the produced offspring is then, E (m) = etp(m, m*)hm n : 

E a (m) = ehrj t ,~ n m"~ 1 m. (CI) 

From iGundersonl |l997l ) we have the yearly mass-specific allocation to reproduction: 

g(m*) = g rjl~ n m"~ 1 , (C.2) 

where go = 1.2 g 1 ~"/year is obtai ned using least sum of squares in fitting the curve to the data for 
oviparous fish in IGundersonl (| 19971 ). Equalling (|C.2|) and E /m allow us to determine the efficiency of 
offspring production e: 

e=^w0.12. (C.3) 



Appendix D. Setting the search rate prefactor 7 from initial feeding level fo 



Food for the smallest individuals in the spectra will be supplied by the background spectrum. If we 
assume that the resource spectrum is at carrying capacity k then an equilibrium initial feeding level fo for 
the small individuals can be calculated using ©. 

Alternatively we may specify an initial feeding level fo and derive one other parameter. By solving 
the feeding level for 7 by using the analytical solution for the density of food <f>(m) (|B.1|) we find 7 as a 
function of fo: 

7 = M „ lnh^l (D1) 

(1 - fo)a^9i, R K (1 - f )y^a6i, R K 



Appendix E. Parameter estimation 

Individual growth: From lKitchell and Stewart] (| 19T7f ) we obtain an estimate of specific dynamic action 
on 15 % of food consumption, and conservative estimates of egestion and excretion on 15 % and 10 % 
respectively. This results in an assimilation efficiency of a = 0.6. 

T he maximum in t ake scales with a 0.6-0.8 exponent, w ith n = 0.75 being an approximate average 
value (jjoblind . 1 19941 ). lAndersen and Riis-Vestergaardl l|2004l ) provides a length-based relationship for the 
maximum int ake rate based on a whiting study adopted for saithe. Using m = O.Oll 3 (m in g and I in cm 
(|Peterd . ll983h , and an energy content of 5.8kJ/g (fish) or 4.2kJ/g (invertebrates) (|Boudreau and Dicki 
Il992l ') yields a prefactor h for the maximal food intake on 83 g 1 ~"/year or 114g 1 - n /year (at 10 ° C). T hese 
intake values overestimate the intake of large individuals since lAndersen and Riis-Vestergaardl (|2004T l use 
an intake exponent of 0.67 instead of n = 0.75. Due to this a value of h — 85g 1_n /year is selected, which 
also provides reasonable fits to 'cod-like' growth curves (m» = 5 kg). 

The standard metabolism scaling expon ent p for fish is slightly higher than for other taxa, around 0.8 
from intra- and interspecies measurements l|Winberel . 1 19561 ; iKillen et all 120071 ) . For simplicity we assume 
p = n. The first term (acquired energy) in the growth model @ is af(m)hm n where the feeding level f(m) 
is a decreasing function of body size (see Results). This has the effect that even when n — p is assumed 
the acquired energy term still effectively scale with a smaller exponent than the maintenance term kmP 
in accordance with the experimental data on food intake and standard metabolism. Furthermore it is 
noted that this clearly makes the individuals in ea ch functional sp e cies n on-neutral. The bioenergetic 
consequences of n 7^ p has been explored in detail by 
The prefactor for standard metabolism can from 



And ersen et all <|2008h . 
Peters! (|l983T Tbe determined to 6.5g 1_n /year if the 



diet is composed of only invertebrates and 4.7 g 1- " /year if all th e energy is from fish. Both values were 
corrected to 10°C using Q10 = 1.83 l|Clarke and Johnston! . 1 1999] ). It is assumed that energy costs due 
to activity can be described with an activity multiplier on the standard metabolic rate. Estimations of 
activity costs are d ifficult to obtain, but activity multi pliers are often reported in the range 1 to 3; e.g. 1.25 
for North Sea cod jHansson et al. . 19961 ). 1.7 for dace jTrudel and Boisclairl. 19961) . and 1 . 44-3.27 for saithe 



Andersen and Riis-Vestergaardl . 



2004) (however see also Rowan and Rasmussenl l| 19961 ); Tile and Steward 



19971 )). A reasonable value of the prefactor for the standard metabolism and activity costs is assumed to 



be k = 10g 1_n /year corresponding to an activity multiplier in the range 1.5 to 2.1. 

Food encounter: The preferred predator-prey mass ratio is set to j3 = 100 (I Jennings et all |2002| ) and 
the width of the selection function to a = 1, which catches the general picture for at least cod a nd dab 
JUrsinl . [l973l ). It should be noted that small organisms such as copepods have a larger a of 3-4.5 (|Ursinl . 
Il974l ). but for simplicity and since focus is on species with rather large m„ the width a will be held 
constant. 

The exponent for swimming speed is q = 0.8 l|Andersen and Beverl . 120061 ). The prefactor 7 for the 
volumetric search rate is difficult to assess from the literature. An alternative approach is to determine it 
as a function of of initial feeding level fo via (|D.1|) . Experience with the model shows that feeding level is 
a decreasing function of body size. This means that it is sensible to use an initial feeding level fo that is 
larger than the expected average feeding level /. In this study a default value of fo = 0.6 is used. This 
along with default parameters yields 7 = 0.8 • 10 4 m 3 g q / yea r (Table U) . An alternative estimate of 7 
is possible by multiplying the prefactors for swimming speed l|Warel . 1 19781 ) and swept reactive field area 
(reactive radius assumed equal to body length): 7 = 20.3 • n ■ 0.01~ 2//3 cm 3 g q /s ~ 4.3 • 10 4 m 3 g 9 /year, 
which indeed justifies the use of fo = 0.6. 

Mortality: Realistic energy reserve sizes may be £ £ [5%; 20%], and in the present study £ = 0.1 is 
used. A widely used background mortality for 'cod-like' m* = 5 kg fishes is Hb = 0.1 year - , which yields 
Mo = 0.84 g 1 "" /year. 



Reproduction: The efficiency of offspring production was not found in the literature. However, an ana- 
lytical expression may be obtained (for n — p) by combining the calculation of yearly mass-sp ecific alloca- 
tion t o reproduction from the bioenergetic model ( |Appendix C| with empirical measuremen ts ([Gunderson 



1997 ). which yields e = Qo/htv 0.1. The fraction of asymptotic size to mature at is n* = 0.25 (lAndersen et al.1 . 
2008). Off spring mass is mo = 0.5 mg which corresponds to an egg diameter of 1 mm (|Woottonl . Il979l ; 
Chamberslll997l ). 



Resource spectrum: The carrying capacity of the resource spectrum sh ould equal the magni t ude o f 
the community spectrum: nm~ x , with an exponent A = 2 — n + q = 2.05 (| Andersen and Bey er. 2006). 
The magnitude of the resourc e spectrum is set to k = 5 • 10 -3 g A_1 /m 3 , which is comparable with 
findings from empirical studies (Rodriguez and Mullin], [l98^) . The constant for resource regeneration rate 
is ro = 4 g 1_p /year at 10° C (I Savage et all [200J). The cut-off of the resource spectrum is set to include 
mesoplankton, m cu t — 0.5 g. 

Appendix F. Expected Lifetime Reproductive Success 

The expected lifetime reproductive success can be split into two components: 1) the probability of 
surviving to become adult, and 2) lifetime reproduction per adult. 



Appendix F.I. Survival probability 

ON 
01 



If we set ^M- = in ([T]) we may obtain the steady-state solution as: 



iV(m) = -^exp(-/^^m], 
g(m,m») \ J g[m,m*) 



where K(m t ) is the constant from the integration along m. We notice that the probability of surviving 
from size mo to size m is p mo -> m = exp ^— g^^m j dm'^ , which along with p mo _> mo = 1 allow us to 
write the survival probability as: 

_ N(m)g(m,m,) 
N(mo)g(m ,m*) 

Appendix F.2. Lifetime adult reproduction 

The amount of energy an adult belonging to a m„ population will spend on reproduction throughout 
its life is: 



Pt„->tV'( m > m*)E(m)dt, 



where t„ is maturation age, and ip(m,mt,)Ei(m) the rate at which energy is allocated to reproduction. 
Noting that g(m,m*) = ^3? allows us to write this as: 

d 1 \ f M ip{m,m*)E{m) 

Riife{m*) = / p m ,^ ra ^- '-^ — -dm. (F.3) 

im, g{m,m*) 



Appendix G. Details of Numerical Methods 

The model is in the form of a series of coupled partial-integro-differential equations (QJ, one for each 
species with the size preference function Q being the integral kernel. The equations are of the first order 
in mass (i.e. hyperbolic) in which case shocks could be formed in the solutions. However the integral kernel 
smooths out any discontinuities and the equations can be solved effectively and accu rately using a standard 
semi- implicit upwind finite-difference scheme for solving PDEs (|Press et all 1199*21 ) . The McKendrick-von 
Foerster PDE |T]) is discretised by calculating g(m,m,) and fi(m) explicitly and making the time update 
implicit, to yield: 

Art+i Mi a i /v i+1 — r/ -, N i+1 

7~7 1 X — —f~w^w , \\ J - i -J 

At Am m 

where i denotes the time step, and w the grid index on the mass axis (i, w £ N + ). First order approximations 
have been used for both the time and mass derivatives. The discretisation in mass is known as the upwind 
approximation since the derivative is calculated from w and w — 1, which is possible since the growth 



function is non-negative. It is further noted that the dm approximation is semi-implicit since densities at 
time step i + 1 are used. Equation (|G.1[) may be written as: 



^:.'\[—p-gi- 1 )+rC 1 (i + -p-gi + Mi^)= v;„. ,g.2) 

1 Aw w J V Aw-w 



which allows us to write an explicit solution for the density spectrum at the i + 1 time step: 

which can be solved iteratively since A^ 1 is given by the boundary condition. The flux in the boundary 
g(mo, m*)N(mo, t) is equal to the reproduction flux of new recruits R (jl If) such that <?o-/Vq +1 = R, which 
yields: A t = 0, and Ci = N{ + -£^R. 

This semi-implicit upwind scheme is very stable but diffusive. The third order QU ICK (Quadrati c 
Upwind Interpolation for Convective Kinematics) scheme along with the techniques by IZiUemal (JT996), 
which prevents overshooting problems introduced by the QUICK method, were used to check that numerical 
diffusion poses no problem for At = 0.02 years, and a m w ma ss grid with 200 l ogarithmically evenly 
distributed points. To ensure stability the Courant condition (i.e. IPress et alj (|l992T l): 

^ < 1, (G.4) 



is prudent to fulfill. The essence of the criterion is that At should be small enough not to allow individuals 
to skip any mass cells m w during their growth trajectory. 

The grid m w spans the offspring size m to 85 kg to include maturation sizes up to the order of 20 kg. 
The grid for the background spectrum ends at m C ut, and the lower limit should be low enough to ensure 
food items for the smallest individuals in the functional species, i.e. 3 decades lower than mo. Identical 
Am m is used for the background and species spectra to ease computations in the overlap [mo; m cut ]- 

To save computational time the ODEs for the background spectrum (If 511 are solved analytically. The 
solution at time to + At for the semi-chemostatic equation is: 

N R (m,t + At) = K(m) - (K{m) - N R (m,t )y- lromF ~ 1+ ^ {m)]At , (G.5) 
where K(m) = is the effective 

carrying capacity at resource size m. 
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